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The heavy quark production in ultraperipheral heavy ion collisions is investigated, with particular 
emphasis on the results from the coherent interactions given by the two-photon process. One 
addresses the heavy quark total cross sections at photon level considering the saturation model 
and the BFKL dynamics in the color dipole picture. The corresponding cross sections at nuclear 
level are presented. It is verified that the QCD dynamics implies an enhancement of the cross section 
in comparison with previous calculations. 
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I. INTRODUCTION 

The behavior of ep/pp scattering in the limit of high center-of-mass energy ^/s and fixed momentum transfer 
is one of the outstanding open questions in the theory of the strong interactions. In the late 1970s, Lipatov and 
collaborators [1] established the papers which form the core of our knowledge of Regge limit (high energy limit) of 
Quantum Chromodynamics (QCD). The physical effect that they describe is often referred to as the QCD Pomeron, 
with the evolution described by the BFKL equation. The simplest process where this equation applies is the high 
energy scattering between two heavy quark- antiquark states, i.e. the onium-onium scattering. For a sufficiently heavy 
onium state, high energy scattering is a perturbative process since the onium radius gives the essential scale at which 
the running coupling a s is evaluated. This process was proposed as a gedanken experiment to investigate the high 
energy regime of QCD in Refs. [2-4] (see also Refs. [5]). In the dipole picture [2], the heavy quark-antiquark pair and 
the soft gluons in the limit of large number of colors N c are viewed as a collection of color dipoles. In this case, the 
cross section can be understood as a product of the number of dipoles in one onium state, the number of dipoles in 
the other onium state and the basic cross section for dipole-dipole scattering due to two-gluon exchange. At leading 
order (LO), the cross section grows rapidly with the energy (a cx a 2 s e( Qjp_1 ) y , where (cvjp — 1) = 4cisNc In 2 and 
Y = In s/Q 2 ) because the number of dipoles in the light cone wave function grows rapidly with the energy. Therefore, 
in principle, the analysis of the energy dependence of the cross section could disentangle the QCD Pomeron effects. 
However, for ep/pp colliders, the study of the QCD Pomeron is made difficult by the fact that the cross section is 
influenced by both short and long distance physics. Only when specific conditions are satisfied is that one can expect 
to determine the QCD Pomeron effects. Some examples are the forward jet production in deeply inelastic events at 
low values of the Bjorken variable x in lepton-hadron scattering and jet production at large rapidity separations in 
hadron-hadron collisions. 

Over the past few years much experimental effort has been devoted towards to observe the QCD Pomeron in 
those processes, but the conclusions are not unambiguous (see e.g. Ref. [6]). This fact is mainly associated with the 
theoretical uncertainty present in the high energy limit of QCD. Only recently, after an effort of ten years, the next- 
to-leading order (NLO) corrections to the BFKL equation were obtained (for a review on NLO BFKL corrections, see 
e.g. Ref. [7] and references therein). Moreover, since the BFKL equation predicts that for s — > oo the corresponding 
cross section rises with a power law of the energy, violating the Froissart bound, new dynamical effects associated 
with the unitarity corrections are expected to stop further growth of the cross sections [8] . This expectation can be 
easily understood: while for large transverse momentum kx, the BFKL equation predicts that the mechanism g — > gg 
populates the transverse space with a large number of small size gluons per unit of rapidity (the transverse size a 
gluon with momentum k± is proportional to l/k'jj, for small k± the produced gluons overlap and fusion processes, 
gg — > g, are equally important. Considering this process, the rise of the gluon distribution with transverse momenta 
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below a typical scale, energy dependent, called saturation scale Q s is reduced, restoring the unitarity. It is important 
to salient that Golec-Biernat and Wusthoff [9] have shown that a saturation model is able to describe the DESY 
ep collider HERA data, in particular the transition from the perturbative to the nonpcrturbative photoproduction 
region, including the inclusive and dim-active cross sections. Therefore, the description of the QCD Pomeron still is 
an open question, which deserves more detailed analysis. 

This situation should be improved in the future with the next generation of linear colliders (For a recent review of 
BFKL searches, see e.g. Ref. [6]). In particular, future e + e _ colliders probably will allow to discriminate between 
BFKL and saturation predictions [10] (For other analysis of photon-photon collisions see Refs. [11-13]). Such reaction 
presents analogies with the process of scattering of two onia discussed above. Although the onium-onium scattering is 
a gedanken experiment, off-shell photon scattering at high energy in e + e~ colliders, where the photons are produced 
from the leptons beams by bremsstrahlung, plays a similar role. In these two-photon reactions, the photon virtualities 
can be made large enough to ensure the applicability of the perturbative methods. Similarly, we can expect to test 
the QCD Pomeron in heavy flavor production in two-photon collisions, where the hard scale is provided by the heavy 
quark mass. 

In the last years many authors have studied in detail the referred process. From the view of the color dipole picture, 
in Ref. [10] the contribution of the dipole-dipole interaction was taken into account considering a generalization of 
the original saturation model [9] to the two-photon reactions. There, the results present good agreement with data 
from e+e~ colliders, remaining some room to discuss issues concerning resolved contributions. In Ref. [14], these 
questions were revisited using the saturation model, where threshold effects were addressed in an accurate way as well 
as the hadronic single-resolved contribution was computed within the dipole formalism. A current open question is 
the discrepancy between the theoretical results and the data on heavy quark production, mostly for the bottom case, 
in the jp photoproduction and two-photon reactions. This feature is also corroborated by the dipole-dipole approach 
from Ref. [15]. As found in Ref. [14], even the inclusion of resolved contributions does not remove completely the 
observed deficit. Quite promising results are obtained in Ref. [16], where heavy quark production in two photon 
reaction is calculated within the fc^-factorization formalism and the unintegrated gluon density for the photon is 
determined using the saturation model. 

^From the view of the BFKL approach, there are several calculations using the leading logarithmic approximation 
[11-13] and considering the next-to-lcading corrections to the total cross section 7*7* process [13,17]. In particular, 
a stable next-to-leading order program relying on the BLM optimal scale setting [19] has produced good results with 
OPAL and L3 data at LEP2 [17]. However, the heavy quark production from real two-photon process is still an open 
question. Only the effects of charm mass in the 7*7* total cross section for the kincmatical region of L3 at LEP was 
analyzed by Bartels et al. [18] considering the LO BFKL dynamics. We address that issue in our analysis, considering 
a BFKL dipole approximation to the QCD dynamics. Similarly to the previous analysis on two-photon physics, one 
shows that the next generation of colliders will be able to discriminate the dynamics associated to the QCD Pomeron. 

Recently, we have proposed in Ref. [20] to investigate QCD Pomeron effects in a different context, namely in photon- 
photon scattering at ultrapcriphcral heavy ion collisions. In this case, the cross sections are enhanced since the 77 
luminosity increases as Z 4 , where Z is the atomic number [21,22]. There, we have analyzed the double diffractive J/^ 
production in 77 collisions, with the photons coming from the Weizsacker - Williams spectrum of the nuclei. For that 
process our results have indicated that future experimental analysis can be useful to discriminate the QCD dynamics 
at high energies. Here we extend our proposal for heavy quark production. In principle, at ultraperipheral heavy 
ion collisions, a photon stemming from the electromagnetic field of one of the two colliding nuclei can interact with 
one photon of the other nucleus (two-photon process) or can penetrate into the other nucleus and interact with its 
hadrons (photon- nucleus process). Therefore, heavy quarks can be produced in photo- nucleus interactions as well as 
photon-photon interactions. The first case has been extensively discussed in literature [23-25] (For a review, see Ref. 
[26]) and recently revised in different approaches [27-29]. On the other hand, heavy quark production in two-photon 
collisions has been discussed essentially in Refs. [29-31]. In these analysis, the underlying contribution comes from 
the QED quark-box diagram (QPM) at LO [Fig. 1 (top panel)], which dominates in the available energies at RHIC, 
W 77 < 6 GeV. In Ref. [29], besides the direct contribution also the resolved contributions were carefully calculated 
and the latter is shown to be negligible. However, at the LHC experiment, where the two-photon energies can reach 
to W 77 < 200 GeV, the QCD diagrams become important and the QCD Pomeron should dominate in this limit [Fig. 
1 (bottom panel)]. Consequently, an enhancement of the cross section associated to the dynamics is expected. We 
will address these issues in the section III, complementing the detailed study presented in Ref. [29]. 

Here, we will restrict our analysis for the two-photon process and its potentiality to investigate the QCD dynamics. 
Similarly as in Ref. [20], our goal here is twofold: to analyze the potentiality of this process to constraint the QCD 
dynamics at high values of energy and to provide reliable estimates for the cross sections concerning that reaction. 
Rclativistic heavy-ion collisions are a potentially prolific source of 77 collisions at high energy colliders. The advantage 
of using heavy ions is that the cross sections varies as Z*a rather then only or. Moreover, the maximum 77 collision 
energy W 77 is 2^/Ra 1 about 6 GeV at RHIC and 200 GeV at LHC, where Ra is the nuclear radius and 7 is the center- 
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of-mass system Lorentz factor of each ion. For two-photon collisions, the cross section for the reaction AA — > AAQQ, 
where Q — c, b, will be given by 



a 



AA^AAQQ 




dT 



'77 
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where t = s/s, s — W 7 is the square of the center of mass (c.m.s.) system energy of the two photons, s of the ion-ion 
system, dC^/dr is the photon luminosity and & qq(s) is the cross section of the 77 interaction (For details related 
to the numerical expressions, see Ref. [20]). Our approach excludes possible final state interactions of the produced 
particles with the colliding particles, allowing reliable calculations of ultrapcripheral heavy ion collisions. Therefore, 
to estimate the heavy quark production it is only necessary to consider a suitable QCD model for photon-photon 
interaction. 

This paper is organized as follows. In Sec. II, we calculate the total cross sections for the process 77 — > QQ based 
on the color dipole picture, using the dipole cross section given by the saturation model and BFKL approach. In Sec. 
Ill, the corresponding heavy quark cross section in coherent ultraperipheral heavy ion collisions are presented, with 
a particular emphasis in the LHC energies. In the last section we drawn our conclusions and summarize the main 
results. 



The heavy quark production via photon-photon collision is a powerful process testing the applicability of perturba- 
tive QCD methods. The photon induced reactions have the advantage of being clear concerning the probing structure 
and the large mass of the heavy quarks produced allows perturbative physics to be employed. The mechanisms of 
production of such process have been systematically studied, despite the available data being modest and having large 
uncertainties (see Ref. [32] and references therein). An important experimental verification is that the data on open 
heavy quark pair production, mostly for bottom quark, stay above theoretical calculations by a sizeable factor. Such 
a feature has generated several studies which rely on higher orders in the collinear factorization approach or on the 
successful results based on the semihard approach. 

An important approach describing two-photon process is given by the color dipole formalism [2,5]. The simple 
physical picture is provided by the photon splitting in a quark-antiquark pair (virtual components in the transverse 
plane), called color dipole, long after the interaction and the further scattering of these pairs with the target. Such an 
approach has produced a consistent and unified description of high energy photon induced processes in deep inelastic 
scattering and two-photon reactions, including the deep inelastic diffractive dissociation. In particular, in the two- 
photon interactions, the target is also given by a color dipole. The wave functions describing the photon virtual 
fluctuation are completely determined from perturbative methods, whereas the effective color dipole cross section has 
to be modelled. It may contain nonperturbative and higher twist contributions as well as parton saturation effects as 
we should discuss latter on. 

The QQ pair production mechanism can be separated into three main contributions: (i) the quark box diagram 
(QPM), with photons coupling to the same heavy quark line, which dominates at not so large photon-photon c.m.s. 
energies, W 77 [Fig. 1 (top panel)]; (ii) the dipolc-dipole interaction, where both photons fluctuate in quark-antiquark 
pairs and interact through gluonic exchange, which dominates at high energies [Fig. 1 (bottom panel)]; (iii) single and 
double resolved contributions, where one of the photons fluctuates into the vector mesons (p, ui, <f>, . . .) interacting 
with the remaining dipole, being a background to the dipole-dipole reaction. The first contribution is well known, 
providing a behavior modulo logarithmic corrections, given by for example on Refs. [33,34]. The second one will 
be reproduced in the following, making use of a sound model for the dipole cross section [10]. The third one will be 
disregarded, but we quote the study presented in Ref. [14], where these pieces were calculated for the first time in the 
color dipole approach. 

The color dipole formulation gives the following expression for the heavy quark pair production in a (real) two- 
photon reactions at high energies, 



where ^ q q } (QQ)( r i z ) are the light (heavy) quark-antiquark wave functions of the photon in the mixed representation. 
The transverse separation (dipole size) of the respective quark pair is denoted by r and their longitudinal momentum 
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fraction denoted by z. The dipole-dipole cross section, cr dd , depends on the dipole sizes and on the parameter x a b 7 
driving the energy behavior. 

First, we write down the expressions for the wave functions squared including the respective quark helicities T, L 
and photon virtuality defined as Q 2 , 

|* t(P| z) \ 2 = Yl 4 ^ + (1 - *f\ 4 K l ( £ / r ) + m ) K o ( £ / ') > ( 3 ) 

/ 

|* i(P| Z )|2 = 4 4 Q 2 , 2 (1 - z) 2 ^ 2 (e f r) , (4) 

/ 

where e/ and m/ stands for the charge and mass of the quark having flavor / and Ko t \ are the McDonald-Bessel 
functions. The auxiliary variable e is defined as e 2 = z(l — z) Q 2 +m'j which in the real photon case, Q 2 = 0, considered 
here, simplifies to s = rrif. Also in such case, the longitudinal photon component is obviously suppressed. The light 
quark masses are considered to be equal and an effective value m q ~ 0.2, is determined from fitting two-photon data 

[10]- 

Lets discuss now a specific model for the dipole-dipole cross section. Here, we follows the formulation of the 
saturation model applied to two-photon reactions [10]. The cross section describing the dipole interaction satisfies 
the saturation property, i.e. a constant dependence on energy at large dipole sizes, and color transparency behavior 
at small dipole sizes. The interface between this two regions is controlled by the saturation scale, which is energy 
dependent and has been determined from fitting data on small- x deep inelastic data [9]. Explicitly, the dipole-dipole 
cross section takes the eikonal-likc form, 



Odd (ri,r 2 ,x ab ) = (T 



1 — cxp 



f 2 



4i? 2 (£„&) 



(5) 



with the overall normalization given by do = (2/3) <r inspired in the quark counting rule and where cr is determined 
for the proton case analysis in Ref . [9] . The effective radius f is constructed is such way to reproduce the dipole-hadron 
parameterization, i.e. f 2 <~ r\ (<~ r 2 .) in the dipole configurations r\ 3> r\ (r\ 3> r 2 ,). Here, we will consider two 
possibilities giving good results as found in Ref. [10]: (I) f 2 — r\r\j(r\ + r 2 ), preferable against two-photon data; 
(II) f 2 = min(r 2 , r^), preferable for charm quark pair production data. The fitted mass for the light quarks takes the 
value m q — 0.21 GeV considering the effective radius (I) and m q = 0.23 GeV using radius (II). The saturation radius 
Ro squared reads as, 

Rl(x ab )= GeV- 2 , (6) 

where the parameters xq and A are also determined in the original saturation model in the proton case [9]. The 
saturation scale, Q s (xab), is obtained by the inverse of the saturation radius described above. The effective variable 
x depends on the two-photon energy and the quark masses and reads as, 

_ 4m 2 +4m 2 

where in the heavy quark production a and b corresponds to a light and a heavy (or vice versa) quark. For low 
energy (large x), the dipole cross section needs to be supplemented by threshold effects. Following [10], they are 
taken into account through the multiplicative factor (1 — x a fc) 5 , coming from dimensional-counting rules. However, a 
different procedure is adopted in Ref. [14], where the threshold is considered by imposing the kinematical constraint 
Mgq + Mqq < W 77 in the integrations at Eq. (2), where M is the invariant mass of the pair. There, also the variable 
x takes a different form than Eq. (7), depending explicitly on z\ and zi- 

In Fig. 2 we reproduce the results on charm production considering the saturation model (SAT-MOD), regarding 
the quark box contribution (QPM) and the dipole-dipole approach using the saturation model. The theoretical 
predictions are compared with the experimental results from the L3 Collaboration [35]. It should be stressed that 
results for the bottom case are straightforward, changing only the quark charge and mass, = 1/3 and = 4.5 
GeV. The charm mass is taken as m c — 1.3 GeV and the two prescriptions for the effective radius are presented. The 
QPM contribution (long-dashed curve) is very important at low energy, starting to be negligible around W 77 = 10 
GeV. Therefore, this contribution should dominate the ultraperiphcral heavy quark production at RHIC energies, 
where the maximal W 77 reached is of order of few GeV's. However, above W 77 > 10 GeV, the dipole-dipole approach 
dominates the physical description, providing an increase of the total cross section. The results using the effective 
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radius (I) (solid line) and (II) (dashed line) are presented. The contribution from dipole-dipole interaction is shown 
separately for the radius (I), represented by the dot-dashed line. This contribution plays a very important role at 
LHC energies, where the two-photon energy in ultraperiphcral collisions can reach up to W ~ 160 GeV. 

Regarding the resolved (single and double) contributions, we have not considered them here. We quote Ref. [14], 
where these contributions are calculated using the saturation model within the dipole-dipole picture. There, it was 
found that for the charm case the single resolved piece constitutes about 30 % of the main contribution. In the bottom 
case, these contributions are quite smaller. 

A prominent challenge in high energy physics concerns the issue what is the correct dynamics driving the interaction 
at asymptotic energies. It is basically dominated by gluonic interaction and the saturation model discussed above 
successfully interpolates between the soft (non-perturbative) and hard (perturbative) regimes, whose transition relies 
on the parton saturation phenomenon. However, the physics of two-photon collisions is an optimal tool testing the 
BFKL approach [1]. In its dipole formulation, the BFKL dynamics is suitable to describe the interaction of small 
size color dipoles from the interacting photons. A lot of work have been done using the LO BFKL approach and 
also attempts to include NLO corrections in two-photon reactions [11-13,17,18]. However, heavy quark production in 
real photon-photon collisions considering the BFKL approach was not addressed in the current literature. Below, we 
present our results for this case. 

To estimate the heavy quark production in the BFKL approach we need a model for the dipole-dipole cross section. 
In particular, the BFKL dipole cross section can be obtained from the momentum representation and it can be 
considered as a solution of the linearized Balitsky-Kovchegov equation in the coordinate space [8,36]. The LO BFKL 
solution is completely known in the momentum space, providing a steep energy dependence associated with the hard 
Pomeron intercept 1 + top and the typical diffusion feature. For our purpose here, we will consider the following 
expression for the dipole-dipole cross section from a BFKL dynamics, 



An 2 ff 2 \ 2 e WQ 3 ln(l/£ ab ) 

3 VI) ^/2inua s ln(l/a; afc ) 



-In 2 (f 2 /f 2 ) 
2ua s ]n(l/xab) 



(I -tab) 5 , (8) 



where a s = N c a s /ir, with number of color iV c = 3. Moreover, for light quarks, C = 1.0, r 2 = 1/Aq CD , where one 
has used A QC d = 0.23. Still, to = 2ip(l) - 2^(1/2) = 4 In 2 and Cj = 28£(3). The Pomeron intercept is given by 
up = 1 + a s u). The prescription used here is not unique. For instance, for dipole-dipole collisions, the choice would 
be f = min(ri, r 2 ) and f = max(n, r 2 ), having impact on the final results. 

In Fig. 3 are presented the results using the dipole-dipole cross section of Eq. (8), considering the effective radius 
(I) and taking into account three different values for the Pomeron intercept. The QPM contribution was also included. 
In our further analysis one considers the intermediate value uip — 1.47, which is steeper than the results from the 
saturation model. It should be stressed that further phenomenology can be done, for instance taking the overall 
normalization, the scale fo and the Pomeron intercept up as free parameters. The differences between the predictions 
of the saturation model and the BFKL approach are already sizeable in the energy interval presented here, with 
important implications in the future TESLA experiment. 

As a final analysis at the photon level, we should calculate the cross section for the interesting final state 2Q 2Q, 
which is by definition an onium-onium scattering. In this case the light quarks in the bottom panel of the Fig. 1 are 
replaced by a pair of heavy quarks. The production of this particular configuration was for the first time addressed 
by Szczurek in Ref. [14] , calling the attention that the dipole-dipole approach gives a unique prediction to production 
of two identical heavy quarks and two identical heavy antiquarks in the final state. The process 77 — ► 2Q 2Q can 
only to be produced in next-to-leading order calculation in the collinear approach and/or hadronization process. The 
cross section for this process is given by, 

a dd ( 77 ^2Q2Q) = J \^ QQ {r u Zl )\ 2 |* Q0 (r 2 , z 2 )| 2 a dd (n, r 2 , x ab ) d 2 r Y d 2 r 2 dz x dz 2 . (9) 

In Fig. 4 are shown the results using the Eq. (9) and considering the saturation model (radius I and II) and the BFKL 
dipole-dipole cross section for the charm quark case. In this case, fg = 1/m 2 and we assume C = 1/9 to reproduce the 
results in the light quark case. Once again, the BFKL result is steeper on energy than the dipole model predictions. 
Similarly to Ref. [14], we obtain that at high energies the cross section for 2c2c is about 9 % of that for the single 
cc production. It is important to salient that in the calculations of the inclusive cc production cross section, the 
contribution of this process should be doubled because each of the quarks can potentially be identified experimentally 
(See discussion in Ref. [14]). In next section, we present our results for heavy quark production in ultraperipheral 
heavy ion collisions, considering both contributions for the inclusive heavy quark production cross section. 
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III. THE HEAVY QUARK PRODUCTION IN ULTRAPERIPHERAL COLLISIONS 



Having determined the cross sections for heavy quark production at the photon level, in the following the results for 
ultrapcriphcral heavy ions collisions will be presented. We give emphasis on the LHC region, where W 77 < 200 GeV, 
and a luminosity C = 4.2 x 10 26 cm _2 s _1 for PbPb collisions (y/s = 5500 GeV). Details on the effective two-photon 
luminosity used here can be found in the previous work in Ref. [20]. As referred above, one has been taken m c = 1.3 
GeV and mb — 4.5 GeV. Further, we discuss the mass dependence of our results. 

Our predictions for the inclusive heavy quark cross section and double cc pair production at \/s = 5500 GeV are 
presented in Table I. In comparison with the predictions from Ref. [29], our results for the inclusive cross section 
are higher by approximately one factor of 1.8. This enhancement is directly associated to the QCD dynamics, which 
resums higher order diagrams beyond those considered previously in the literature. It is important to salient that 
this factor is larger than the uncertainty associated with the treatment of the 77 luminosity, which is estimated to 
be about 15% [29]. For RHIC energies (not presented in the Table), due to the limited center of mass energy of the 
photon-photon process, the QCD Pomeron effects are not important, and our results are similar with those presented 
in Ref. [29]. 

In Fig. 5 are presented the results for the energy dependence of the inclusive charm pair production in two-photon 
ultraperipheral heavy ion collisions. The estimates using the saturation model are shown as the solid and long-dashed 
lines, corresponding to the choice (I) and (II) for the effective radius, respectively. The dipole BFKL calculation 
is shown in the dot-dashed curve, where we have considered the Pomeron intercept ujp = 1.47. For comparison, 
the QPM and the single cc cross sections are also presented. Despite the deviations between the approaches being 
large at the photon level, the results for the peripheral reaction are quite similar. This feature is expected, since 
the two-photon luminosity in nuclear collisions strongly suppress high energy contributions in comparison with the 
behavior present in e + e~ collisions. Concerning the quark mass dependence of our results, we have checked that if 
one uses m c = 1.2 GeV, as done in [29], the cross section grows approximately 25 %. Therefore, for this quark mass 
choice, we predict an enhancement of a factor two. Here we have used m c = 1.3 GeV, since with this choice the 
experimental data are well described in the dipole approach. 

In Fig. 6 are shown the results for the energy dependence of the inclusive bottom pair production in two-photon 
ultraperipheral heavy ion collisions. The solid line represents the estimate obtained with the saturation approach 
(model I), the dot-dashed line the dipole BFKL result and the dashed line is the QPM prediction. The deviation 
between the saturation and BFKL approaches reaches up to 10%. Similarly to the charm, the analysis of the bottom 
quark production in ultraperipheral heavy ion collisions would not allow to discriminate between the distinct kinds of 
dynamics. However, these effects cannot be disregarded at LHC energies, since our predictions are a factor two larger 
than those obtained in [29]. 

A comment is in order here. In [29] a detailed analysis of the experimental separation between photoproduction 
and two-photon interactions was presented. There, the authors have estimated that the two-photon cross sections are 
at least 1000 smaller than the photoproduction cross section, which makes the experimental separation between the 
two interactions very hard. Our calculations indicate that the inclusion of the QCD Pomeron effects implies higher 
cross sections at two-photon level and, consequently, larger cross sections in ultrapcriphcral collisions. Therefore, 
the inclusion of these effects implies that the contribution of heavy quark production in two-photon interactions is 
non-negligible. However, the experimental separation of two-photon process is still remains a challenge. 

In Fig. 7 we present separately the results for the double charm pair production in two-photon ultraperipheral 
heavy ion collisions, obtained from Eq. (9). The solid line represents the estimate obtained with the saturation 
approach (model I) and the dot-dashed line is the dipole BFKL result. The difference between the predictions of the 
two approaches is about 20% in this particular case. In comparison with inclusive charm pair production, the cross 
section is two order of magnitude smaller, but still sizeable for an experimental analysis. We salient that the efficiency 
of flavor tagging is in general very small, which also makes the experimental detection of this process a challenge. 
However, it important to emphasize that distinctly from the single pair production which is one background of the 
photoproduction process, the double heavy quark production can be produced in photo-nucleus reactions only if the 
multiple interactions are considered. Since, in general, these interactions are suppressed by a factor of about 1000 
[37], we have that for this process two photon interactions and photon- nucleus interactions would be of same order. 

IV. SUMMARY 

We have investigated the heavy quark production in the ultraperipheral heavy ion collisions from two-photon 
reactions. The total cross sections at the photon level are calculated using the saturation model, describing the 
dipolc-dipole interaction at high energies including in a suitable way the saturation phenomenon. The dipole BFKL 
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approach is also considered as an important test of QCD dynamics and predictions for the photon-photon cross section 
using this approach are presented. Both calculations are in reasonable agreement with the current experimental 
measurements. In particular, we have considered a typical final state allowed by the color dipole picture, namely the 
double heavy quark pair production 77 — * 2Q 2Q. 

The resulting cross sections in ultrapcripheral heavy ion collisions are estimated in these distinct dynamical ap- 
proaches for the QCD Pomeron and presented. Our results demonstrate that the analysis of the heavy quark produc- 
tion in ultraperipheral heavy ion collisions would not allow to discriminate between the distinct dynamics. However, 
the QCD Pomeron effects cannot be disregarded at LHC energies, since the cross sections are enhanced by a factor 
two for these effects. Moreover, we have analyzed the double charm pair production and verified that the cross sec- 
tion for this process is sizeable, mainly due to the fact that photoproduction background can only be produced by 
higher-order multiple interactions. Our results motivate more detailed studies. In particular, the analysis of the mass 
and transverse momentum distributions could be useful to discriminate between photon-Pomeron and two-photon 
processes. Work in this direction is in progress. 
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SAT-MOD (I) 


SAT-MOD (II) 


LO BFKL 


cc 


1810000 nb 


1900000 nb 


1951200 nb 


bb 


2000 nb 




2200 nb 


2c 2c 


40000 nb 




48000 nb 



TABLE I. The heavy quark production total cross sections for ultraperipheral heavy ion collisions at LHC (y/s = 5500 GeV) 
for PbPb. 
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FIG. 1. The diagrams contributing to the QPM piece (top panel) and dipole-dipole interactions (bottom panel). Figures 
from Ref. [10]. 
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FIG. 2. The cross section for the inclusive charm production in two-photon reactions. The QPM contribution (thin dashed 
line) is shown separately. The saturation model considering different effective radius are represented by the solid (model I) and 
long-dashed (model II) curves. Data from L3 Collaboration. 
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FIG. 3. The cross section for the inclusive charm production in two-photon process considering the BFKL dipole-dipole 
approach. The results for three different values for the Pomeron intercept top are presented. The QPM contribution is also 
included. Data from L3 Collaboration. 
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FIG. 4. The cross section 17(77 — * 2c 2c) considering the saturation model (models I and II) and the BFKL approach, with 
a Pomeron intercept ujp — 1.47. 
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FIG. 5. Inclusive charm pair production in ultraperipheral heavy ion collisions considering the distinct QCD approaches for 
the 77 interaction. For comparison the QPM (dotted line) and single cc (dashed line) cross section are also presented. In the 
last case, the contribution of double cc production for the inclusive cross section was subtracted. 
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FIG. 6. Inclusive bottom pair production in ultraperipheral heavy ion collisions. For comparison the QPM prediction is also 
presented. 
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FIG. 7. Double pair charm production in (two-photon) Ultraperipheral Heavy Ion Collisions. 
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